Circular motion of asymmetric self-propelling particles 
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Micron-sized self-propelled (active) particles can be considered as model systems for characterizing 
more complex biological organisms like swimming bacteria or motile cells. We produce asymmetric 
microswimmers by soft lithography and study their circular motion on a substrate and near channel 
boundaries. Our experimental observations are in full agreement with a theory of Brownian dynamics 
for asymmetric self-propelled particles, which couples their translational and orientational motion. 
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Micron-sized particles undergoing active Brownian mo- 
tion [1 currently receive considerable attention from ex- 
perimentalists and theoreticians because their locomo- 
tion behavior resembles the trajectories of motile mi- 
croorganisms [2^. Therefore, such systems allow in- 
teresting insights into how active matter [6] organizes 
into complex dynamical structures. During the last 
decade, different experimental realizations of microswim- 
mers have been investigated, where, e.g., artificial flag- 
ella or thermophoretic [8] and diffusiophoretic [9] 
driving forces lead to active motion of micron-sized ob- 
jects. So far, most studies have concentrated on spher- 
ical or rod-like microswimmers whose dynamics is well 
described by a persistent random walk with a transition 
from a short-time ballistic to a long-time diffusive behav- 
ior [To]. Such simple rotationally symmetric shapes, how- 
ever, usually provide only a crude approximation for self- 
propelling microorganisms, which are often asymmetric 
around their propulsion axis. Then, generically, a torque 
is induced that significantly perturbs the swimming path 
and results in a characteristic circular motion. 

In this Letter, we experimentally and theoretically 
study the motion of asymmetric self-propelled particles 
in a viscous liquid. We observe a pronounced circular mo- 
tion whose curvature radius is independent of the propul- 
sion strength but only depends on the shape of the swim- 
mer. Based on the shape-dependent particle mobility 
matrix, we propose two coupled Langevin equations for 
the translational and rotational motion of the particles 
under an intrinsic force, which dictates the swimming ve- 
locity. The anisotropic particle shape then generates an 
additional velocity-dependent torque, in agreement with 
our measurements. Furthermore, we also investigate the 
motion of asymmetric particles in lateral confinement. 
In agreement with theoretical predictions we find either 
a stable sliding along the wall or a reflection, depending 



on the contact angle. 

Asymmetric L-shaped swimmers with arm lengths of 
9 and 6 |Lim were fabricated from photoresist SU-8 by 
photolithography [11]. In short, a 2.5 |LLm thick layer 
of SU-8 is spin coated onto a silicon wafer, soft-baked 
for 80 s at 95 °C and then exposed to ultraviolet light 
through a photo mask. After a post-exposure bake at 
95 °C for 140 s the entire wafer with the attached par- 
ticles is coated with a 20 nm thick Au layer by thermal 
evaporation. When the wafer is tilted to approximately 
90° relative to the evaporation source, the Au is selec- 
tively deposited at the front side of the short arms as 
schematically shown in Figs.[TJa),(b). Finally, the coated 
particles are released from the wafer by an ultrasonic 
bath treatment. A small amount of L-shaped particles is 
suspended in a homogeneous mixture of water and 2,6- 
lutidine at critical concentration (28.6 mass percent of 
lutidine), which is kept several degrees below its lower 
critical point (Tc = 34.1 °C) [12 . To confine the par- 
ticle's motion to two dimensions, the suspension is con- 
tained in a sealed sample cell with 7 iitm height. The 
particles are localized above the lower wall at an average 
height of about 100 nm due to the presence of electro- 
static and gravitational forces. Under these conditions, 
they cannot rotate between the two configurations shown 
in Figs. [TJa),(b), which will be denoted as L+ (left) and 
L- (right) in the following. When the sample cell is illu- 
minated by light (A = 532 nm) with intensities ranging 
on the order of several M-W/ jitm^, the metal cap becomes 
slightly heated above the critical point and thus induces 
a local demixing of the solvent [131 Ej- This leads to 
a self-phoretic particle motion similar to what has been 
observed in other systems [ISWTT]. 

Figures [TJa),(b) show trajectories of L+ and L- swim- 
mers obtained by digital video microscopy for an illu- 
mination intensity of 7.5 jLtW/iLtm^, which corresponds 
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Figure 1. (Color online) (a),(b) Trajectories of an (a) L+ and 
(b) L- swimmer for an illumination intensity of 7.5 |j,W/^Lm^. 
(Red) bullets and (blue) square symbols correspond to ini- 
tial particle positions and those after 1 min each, respectively 
The insets show microscope images of two different swimmers 
with the Au coating (not visible in the brightfield image) indi- 
cated by an arrow. (c),(d),(e) Probability distributions p(a) 
of the angle a (see inset in (c)) between the normal vector 
of the metal coating and the displacement vector Ar of an L+ 
particle in time intervals of 12 s each for illumination intensi- 
ties (c) I = ^lW/^lm^ (d) 5 ^lW/^lm^ and (e) 7.5 [iW/[im^. 



to a mean propulsion speed of 1.25 |Lim/s. In contrast 
to spherical swimmers, here a pronounced circular mo- 
tion with clockwise (L+) and counter-clockwise (L-) di- 
rection of rotation is observed. For the characterization 
of trajectories we determined the center-of-mass position 
r(t) = {x{t)^y{t)) and the normalized orientation vector 
u_L of the particles (see inset of Fig.jljc)). From this, we 
derived the angle a between the displacement vector Ar 
and the particle's body orientation u^. Figures [TJc)-(e) 
show how the normalized probability distribution p{a) 
changes with increasing illumination intensity /. In case 
of pure Brownian motion (see Fig. [ijc)) p{a) ^ const, 
since the orientational and translational degrees of free- 
dom are decoupled when only random forces are acting 
on the particle. In presence of a propulsion force which 
is constant in the body frame of the particle, however, 
the translational and rotational motion of an asymmet- 
ric particle are coupled. This leads to a peaked behavior 
of p{a) as shown in Figs. [TJd),(e). The peak's halfwidth 
decreases with increasing illumination intensity since the 
contribution of the Brownian motion is more and more 
dominated by the propulsive part. In addition, the peaks 
are shifted to positive (negative) values for a particle 
swimming in (counter-) clockwise direction. The posi- 
tion of the peak is given by a = 7rAt/r, where r is the 
intensity-dependent cycle duration of the circle swimmer 
(cf.. Fig. |2|b)) and At is the considered time interval. 
This estimate (see arrows in Figs. [TJd) (r = 60 s) and 
[TJe) (r = 40 s)) is in good agreement with the experi- 
mental data. The shift of the maximum of p{a) docu- 
ments a torque responsible for the observed circular mo- 
tion of such asymmetric swimmers. In contrast to an 



externally applied constant torque [18], here it is due to 
viscous forces acting on the self-propelling particle. This 
is supported by the experimental observation that the 
particle's angular velocity uj{t) = da/dt increases lin- 
early with its total translational velocity v{t) (see Fig. 
|2|a)). As a direct result of the linear relationship be- 
tween u and V, the radius R of the circular trajectories 
becomes independent of the propulsion speed, which is 
set by the illumination intensity (see Fig. [2|b)). 

For a theoretical description of the motion of asymmet- 
ric swimmers, we consider an effective propulsion force F 
[19], which is constant in a body- fixed coordinate sys- 
tem that rotates with the active particle. With the unit 
vectors u_\_ = (— sin^, cos^) and uy = (cos 0, sin 0) (see 
Figs. [TJc) and [3]^ a)), where - in case of L-shaped parti- 
cles - (j) is the angle between the short arm and the x 
axis, the propulsion force F is given by F = Fuint with 
Uint = (cu|| + ux)/\/l + with the constant c depend- 
ing on how the force is aligned relative to the particle 
shape. If the propulsion force is aligned along the long 
axis u^, one obtains c = 0, i.e., uint = u±. In case of 
an asymmetric particle, the propulsion force leads also 
to a velocity-dependent torque relative to the particle's 
center-of-mass. For c = this torque is given hy M = IF 
with / the effective lever arm (see Fig. |3|a)). Our the- 
oretical model is valid for arbitrary particle shapes and 
values of c and However, for the sake of clarity, we 
set c = as this applies for the L-shaped particles con- 
sidered here. Accordingly, we obtain the following cou- 
pled Langevin equations, which describe the motion of 
an asymmetric microswimmer 



r = /3F(Dtu± + /Dc) 
(/) = /3F(/I)r + Dc-u±) 



(1) 



Here, P = 1/(/cb^) is the inverse effective thermal en- 
ergy of the system. These Langevin equations contain 
the translational short-time diffusion tensor Dt(0) = 
I)||U|| (g)U|| +I)^(u|| (g)ux + ux0U||) + l^±u±0ux with the 
dyadic product (g) and the translation-rotation coupling 
vector Dc(^) = D^c^\\ ^ D^u±_ [20]. The translational 
diffusion coefficients Z^y, D^^ and D±^ the coupling coef- 
ficients and and the rotational diffusion constant 
are determined by the specific shape of the particle. 
Finally, Cr(^) C^l^) Gaussian noise terms of zero 
mean and variances (Cr(^i) ^ Cr(^2)) = 20x^(^1 — ^2), 
(Cr(^i)C0fe)) = 2Dc^(ti -^2), and ((^(tOC^fe)) = 

In case of vanishing noise, Eq. ([T]) immediately leads 
to a circular trajectory with radius 



R = 



\ 



{D^ + ID^Y 



(2) 



In agreement with the experimental observation (see 
Fig. |2|b)) the radius does not depend on the parti- 
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Figure 2. (Color online) (a) Angular velocity cj and (b) ra- 
dius R of the circular motion of an L+ swimmer plotted as 
functions of the linear velocity = |v| and the illumination 
intensity I ^ v. The dashed lines correspond to a linear fit 
with nonzero and zero slope, respectively. 



cle velocity set by the propulsion force. Rather, the 
value of R is only determined by the particle's geom- 
etry, which defines its diffusional properties. More- 
over, the translational and angular particle velocities 

are v = pF^{Df + -j- {D^ + ID^f and uj = 

PF{Dq + ID^i). Both quantities are proportional to the 
internal force F and ensure R = v/\uj\ in perfect agree- 
ment with the experimental results shown in Fig. [2|a). 

For a quantitative comparison with the experimental 
data, most importantly, the diffusion and coupling co- 
efficients for the particles under study have to be de- 
termined. They constitute the components of the gen- 
eralized diffusion matrix and are, in principle, obtained 
from solving the Stokes equation that describes the low 
Reynolds number flow field around a particle close to the 
substrate [22^. This procedure can be approximated by 
using a bead model where the L-shaped particle is 
assembled from a large number of rigidly connected small 
spheres. Exploiting the linearity of the Stokes equation, 
the hydrodynamic interactions between any pair of those 
beads can be superimposed to calculate the generalized 
mobility tensor of the L-shaped particle and from that 
its diffusion and coupling coefficients; details of the cal- 
culation are outlined in Ref. [23] . This method is well es- 
tablished for arbitrarily shaped particles in bulk solution 
[23, 24 . We take into account the presence of the sub- 
strate by using the Stokeslet close to a no-slip boundary 
[25j to model the hydrodynamic interactions between the 
component beads in the bead model. For the L-shaped 
particles considered here, we find that the value of D_i 
exceeds the terms Dj-^ ^^c^ numerator 
of Eq. (H by more than one order of magnitude (given 
that / is in the range of 1 |LLm). On the other hand, the 
value of is negligible compared to IDr. This finally 
yields 



R=\D^/{lDn)\ 



(3) 



as an approximate expression for the trajectory radius 
and, correspondingly. 



for the angular velocity. 



Table L Diffusion coefficients for the L-shaped particle in Fig. 
^a) on a substrate: translational diffusion along the long 
\D±) and the short {D\\) axis of the L-shaped particle as well 
as rotational diffusion constant Dr. 
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We determined the diffusion coefficients D^^ D\\^ and 
I^R experimentally under equilibrium conditions (i.e., in 
the absence of propulsion) from the short-time correla- 
tions of the translational and orientational components of 
the particle's trajectories [26','27 (see Tab.|l|. The exper- 
imental values are in good agreement with the theoretical 
predictions. 

Inserting the experimentally determined values for 
the diffusion coefficients and the mean trajectory ra- 
dius R = 7.91 |Ltm into Eq. (|3|, we obtain the effective 
lever arm / = —1.65 iitm. This value is about a factor 
of two larger compared to an ideally shaped L-particle 
(see Fig.jsja)) with its propulsion force perfectly centered 
at the middle of the Au layer. This deviation suggests 
that the force is shifted by 0.94 jam in lateral direction, 
which is most likely caused by small inhomogeneities of 
the Au layer due to shadowing effects during the graz- 
ing incidence metal evaporation. Accordingly, from Eq. 
Q we obtain the intensity-dependent propulsion force 
F/I = 0.83 X 10-^^ N^imV^tW. 

To compare the trajectories obtained from the 
Langevin equations ([T]) with experimental data, we di- 
vided the measured trajectories into smaller segments 
and superimposed them such that the initial slopes and 
positions of the segments overlap. After averaging the 
data we obtained the noise-averaged mean swimming 
path, which is predicted to be a logarithmic spiral {spira 
mirabilis) [28 that is given in polar coordinates by 



r(0) = pF. 



■ exp 



— ( 

UJ 



(5) 



Lo = /3DrIF 



(4) 



Qualitatively, such spirals can be understood as follows: 
in the absence of thermal noise, the average swimming 
path corresponds to a circle with radius R given by 
Eq. (|3|. In the presence of thermal noise, however, 
single trajectory segments become increasingly different 
as time proceeds. This leads to decreasing distances 
di between adjacent turns of the mean swimming path 
{di/di-^i = exp {27rDji/\uj\), see Fig.jsj^d)) and, finally, to 
the convergence in a single point for t ^ oo. Due to the 
alignment of the initial slope, this point is shifted relative 
to the starting point depending on the alignment angle 
and the circulation direction of the particle. 
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Figure 3. (Color online) (a) Geometrical sketch of an ideal 
L+ swimmer as considered in our model. The dimensions are 
a = 9 ^im, 6 = 6 ^im, xs — 2.29 p.m, and ys — 3.55 p.m (for 
homogeneous mass density and an additional 20 nm thick Au 
layer) . The internal force F induces a torque M on the center- 
of-mass S depending on the lever arm /. (b),(c) Visualization 
of the experimental trajectory (for an illumination intensity of 
/ = 7.5 M^W/ \iv£?) that is used for the quantitative analysis of 
the fluctuation-averaged trajectory in (d). The dashed curve 
in (d) is the experimental one, and the solid curve shows the 
theoretical prediction with the starting point indicated by a 
red bullet. Inset: close-up of the framed area in the plot. 



The solid curve in Fig. [sj^d) is the theoretical predic- 
tion (see Eq. ([5|) with the measured values of D I)r, 
and uj. On the other hand, the dashed curve in Fig. 
[sj^d) visualizes the noise- averaged trajectory determined 
directly from the experimental data (see Figs. [3]^b),(c)). 
The agreement of the two curves constitutes a direct ver- 
ification of our theoretical model on a fundamental level. 

Finally, we also address the motion of asymmetric 
swimmers under confinement, e.g., their interaction with 
a straight wall. This is shown in Fig. [4]^ a) exemplarily for 
an L+ swimmer which approaches the wall at an angle 
d. Due to the internal torque associated with the active 
particle motion, it becomes stabilized at the wall and 
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Figure 4. (Color online) (a),(b) Trajectories of an L+ swim- 
mer approaching a straight wall at different angles (symbols 
correspond to positions after 1 min each), (c) Experimen- 
tally determined particle motion for different contact angles 
Q. Bullets and open squares correspond to particle sliding and 
reflection. (d),(e) Visualization of the predicted types of mo- 
tion for an L+ swimmer with arrows indicating the direction 
of the propulsion force: (d) stable sliding and (e) reflection. 
The angles are defined in the text. 



smoothly glides to the right along the interface. In con- 
trast, for a much larger initial contact angle the internal 
torque rotates the front part of the particle away from 
the obstacle, the motion is unstable, and the swimmer 
is reflected by the wall (see Fig. |4jb)) [29 . Figure [4]^ c) 
shows the observed dependence of the motional behavior 
as a function of the approaching angle. 

The experimental findings are in line with an instabil- 
ity analysis based on a torque balance condition of an 
L-shaped particle at wall contact as a function of its con- 
tact angle d. For 6>crit < 6> < tt (see Figs. |4jb),(e)) with 
a critical angle ^crit, the particle is reflected, while for 
^crit (see Figs. |4]^a),(d)) stable sliding with an 
angle ^gi occurs. Both, and 6>crit are given as stable 
and unstable solutions, respectively, of the torque bal- 
ance condition 



\(a — ys) cos — xs sin 6] sin . 



(6) 



For / = —0.71 |Ltm corresponding to an ideal L-shaped 
particle with the propulsion force centered in the middle 
of the Au layer, we obtain ^gi = 8.0° and ^cnt = 59.2°, 
which is in good agreement with the measured value of 
about ^crit = 60° (see Fig. [4]^c)). The observed scatter in 
the experimental data around the critical angle is due to 
thermal fluctuations that wash out the sharp transition 
between the sliding and the reflection regime. 

In conclusion, we have demonstrated that due to 
viscous forces of the surrounding liquid, asymmetric 
microswimmers are subjected to a velocity-dependent 
torque. This leads to a circular motion, which is ob- 
served in experiments in agreement with a theoretical 
model based on two coupled Langevin equations. In a 
channel geometry, this torque leads either to a reflection 
or a stable sliding motion along the wall. An interesting 
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question for the future addresses how asymmetric swim- 
mers move through patterned media. In the presence of 
a drift force, one may expect Shapiro steps in the particle 
current similar to what has also been found in colloidal 
systems driven by a circular drive [30^. Another appeal- 
ing outlook addresses the motion of chiral swimmers in 
the presence of external fields such as gravity [31 . In 
case of asymmetric particles, this leads to an orienta- 
tional alignment during their sedimentation, which may 
result in a preferential motion relative to gravity similar 
to the gravitactic behavior of asymmetric cells as, e.g., 
Chlamydomonas [32l [33]. 
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